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A one parameter solvable model for three bosons subject to delta function attractive interactions 
in one-dimension with periodic boundary conditions is studied. The energy levels and wave func- 
tions are classified and given explicitly in terms of three momenta. In particular, eigenstates and 
eigenvalues are described as functions of the model parameter, c. Some of the states are given in 
terms of complex momenta and represent dimer or trimer configurations for large negative c. The 
asymptotic behaviour for small and large values of the parameter, and at thresholds between real 
and complex momenta is provided. The properties of the potential energy are also discussed. 



I. INTRODUCTION 

Quantum solvable models where the wave functions, energy eigenvalues and other quantities of physical interest 
can be obtained explicitly in terms of known functions or with minimum numerical effort (typically by solving a 
transcendental equation or by quadratures) are useful to test and refine concepts and/or numerical methods, and as 
first approximations to more realistic systems. Occasionally unexpected physical phenomena are revealed 0. In this 
paper we shall analyze a one parameter model for three bosons subject to attractive delta function pair interactions 
in one-dimension with periodic boundary conditions (contrast this to three particles "on a ring", see Q). Our original 
motivation was to examine a system with attractive forces where single and compound particles may coexist. This is 
of particular interest when studying the kinetic theory of gases composed by particles that can form stable aggregates 
(such as dimers or trimers), especially when chemical reactions can occur p]. This article deals exclusively with the 
model itself, which has been found to be quite complex in mathematical detail. 

The literature on one-dimensional solvable models of three, and generally N, particles is rather extensive. These 
models could be primarily classified according to the type of interaction involved However, even with the same 
interaction but with different boundary conditions, different formal treatments are required and very different results 
may be found. Periodic boundary conditions are suitable for modelling a gas or a crystal lattice in the thermodynamic 
limit. In contrast, in the limit as the box length becomes large, information about the corresponding scattering 
problem of a one-dimensional system can be extracted ^]J^). Indeed these are the standard boundary conditions 
for actual calculations of time dependent wave function scattering [Q. The model studied here is a particular case 
of the "interacting bose gas" of Lieb and Liniger ||, who examined N particles subject to 2-body delta function 
interactions and boson symmetry. Further analysis of this gas was carried out in several papers H [L2), but, having 
different objectives and applications in mind, in none of these works was the attractive case examined, except for the 
appendix on the N — 2 case in || . Lieb and Liniger found some unexpected effects of the periodic confinement but 
did not investigate the analogous effects for N > 2. In a series of papers ]f^I5t, McGuire has examined a related one- 
dimensional many-particle Fermion system with one particle having spin down in a sea of spin up particles, interacting 
via (5-function potentials (both repulsive and attractive). For other models with periodic boundary conditions but 
different interactions, see []1 6| — |20f1 . 

If the particles are not confined in a box, the wave function obeys the standard vanishing boundary conditions 
(for bound states) or scattering boundary conditions, at infinite distances However, in most available models, 
rearrangement processes where a bound pair collides with a single particle to form a new pair are not allowed and 
cannot be examined. An exception is the work of McGuire on the attractive, two body, delta function interaction 
describing the scattering wave functions and the bound (A^-body) states |22|. This model has been generalized, 
examined by means of several formalisms, or applied for different purposes p3[-p8|]. As stated before, the attractive 
case for bosons has not been examined with periodic boundary conditions and the present work fills this lacuna for 
N = 3. 
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II. MODEL DESCRIPTION 



The stationary Schrodinger equation for three equal mass particles in 1-dimension with 2-body delta function 
interactions reads 

-£eS+ 25 E%*-^=^, w 

i=l y * i<j 

where y t (i = 1, 2, 3) are the particle coordinates. If they are enclosed in a box with length L it is convenient to divide 
this equation by h 2 /(2mL 2 ) and use instead 

3 Q^Jj 

-Et^ + ^E^-^ = (2) 

expressed in terms of the dimensionless quantities 

Xi = yi/L (3) 
E = 2mL 2 E/h 2 (4) 
c = 2mcL/h 2 . (5) 

Note that c is the only parameter of the model, and that, according to (|J), the effect of enlargening the box or making 
the interaction stronger are equivalent. This work addresses the attractive case in particular, corresponding to c < 0. 
But the repulsive case (c > 0) is also covered since the nature of the solutions for the repulsive case are the same 
as one class of solutions for the attractive case. Of particular interest is how the wavenumbers associated with a 
wavefunction change, as the potential parameter c varies continuously from repulsion to attraction. 

The delta function potential produces a jump in the derivatives of the wave function where two particles meet. 
This jump is proportional to c and to the wave function at that point, 

= 2cip\ Xj = Xk . (6) 

Since the delta function interaction allows the particles to cross each other, all orderings are possible, each ordering 
corresponding to a particular "region" of coordinate space and one of the permutations of the three particles. On the 
basis that the particles are bosons, it is sufficient to study the wave function in only one of these regions, specifically 
the "primary" region 

#123 : < Xi < x 2 < x 3 < 1. (7) 

The wave function in any other region Rijk of coordinate space is then simply obtained from the wave function in 
R123 by interchanging the particle labels. In region -R123, the equations (||) and (||) can be written as 

3 Q2 







(dip 


_ _c¥ \ 


\dxj 


9xk) Xj=Xk+ 


\dxj 


dx k ) 



E ^ = Ei> (8) 



; = 1 

for xi 5^ a; 2 7^ X3 =/= x±, and the jump conditions 

3 -^1* 



dx.j + i dxj 



— cip 



(9) 

Xj + i—Xj 



If periodic boundary conditions are also imposed, a displacement of 1 in any of the coordinates Xj leaves the function 
unchanged. In particular, in -R123 this means that the wave function obeys 

1p(0,X2,X 3 ) = ip(x 2 ,x 3 ,l). (10) 

Similarly, the derivatives satisfy 

d d 

■7^1p(x,X2,X 3 )\ x =o = -Q^^{X2,X 3 ,X)\ X = 1 . (11) 
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In a two-body collision between particles of equal mass, the delta function interaction can only interchange the 
momenta of the incident particles or leave them unchanged. In other words, there is no diffraction, that is, no "new" 
momenta, different from the initial ones, are created. For three equal mass particles with delta two-body interactions, 
the eigenstates can thus be written in terms of only three plane waves with (dimcnsionless) "momenta" {kj}, 



ip(xi < X2 < ^3) = a(P)P exp 
p 




(12) 



where the sum is over all permutations P of the {kj}, and a(P) are coefficients to be fixed by the boundary conditions 
determined by the periodicity, Eqs. ( |lCi| ) and (|Tl|), and by the delta function interaction, Eq. (||). This wave function 
form is known as the "Bethe ansatz" and it was first applied to spin chains [^9) . In the context of the Bose gas, Yang 
and Yang []l2| used a continuity argument to show that for positive c all states are given by ([l2]) with real fc's, a result 
which was later established rigorously by Dorlas [30| . In the present work it is shown that all states for c > are 
continuously connected to c < states in k space so we are confident that the eigenstates discussed later form in fact 
a complete set. An important difference with the repulsive case is that for c < the fc's may become complex. 
From (|J) and (|lj) the energy is simply obtained as 



E : 



(13) 



but this should not be interpreted as purely kinetic energy since there is generally a potential energy contribution 
to E. Note that Eq. (||) is the Schrodinger equation only when the positions of the particles are all different. The 
true kinetic energy has to take into account the jumps in the wave function derivative at the region boundaries. The 
calculation of the potential energy is however somewhat involved and is discussed in the Appendix. 

The structure of the coefficients a(P) is imposed by the jump boundary condition (|^) as explained, e.g., in The 
amplitudes for two permutations differing by a transposition of two particles are related by a factor —e lSje , 

a(123) = 1 
a(213) = ~e lf>21 
a(132) = -e 1032 
a(321) = _ e i(e21+031+e32) 
a(312) = 

a(231) = e ^+ e 3i) 



where 



i(kj - ki) 
c + i(kj — kf) 



(14) 
(15) 



(16) 



By substituting ( |l2|) and ( |l4| ) into the periodicity conditions (|l0| ) and (|Tl]), the following set of coupled transcendental 
equations is found 



exp 



.7 = 1,2,3. 



(17) 



where, by convention, 9jj = 0. Solving Eq. (h_7j) for the kj gives 



2nmi 
27rm 2 
27rm 3 



721 
?12 
?23 



^3 1 
^32 
?13, 



for some set of integers {nij}, while the 9ji are given in terms of the kj by 

c + i(kj — kf) 



i In 



i(kj 



= — 2 arctan 1 



(18) 



(19) 
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While it may appear to be natural to choose the principal branch of the logarithm and the arctangent, another choice 
is more appropriate. Since a state is uniquely defined by the set of numbers {kj}, irrespective of the order, because 
the particles are bosons, it is convenient to order the {kj}, when they are real, according to 



k\ < k 2 < k 3 . 

Consistent with this, the ranges of 621, #32 and #31 are chosen to satisfy 



2tt < 



W < 0, 



(20) 



(21) 



on the basis that the corresponding kj — kg are positive. In this way, solving for the set of {kj} is equivalent to solving 
for the set of dji and it is noted that, for the above choice for the range of 9j£ , the 9j£ vary continuously as c and/or 
the fcj's vary continuously over their allowed ranges. This is also true if some of the kj and/or 9jt become complex, 
as discussed in latter sections. Thus the rrij provide a unique classification of the energy levels, and for a given set 
of rrij , the energy eigenvalue (and eigenvector) can be followed continuously as a function of c, as c varies from 00 to 
—00. 

On taking the product of the three equations of the form of Eq. (|l7|), it follows that the total momentum p is 
quantized, 



P 



kj = 2i\(m\ + 777,2 + Tris) = 2im p 



(22) 



This is the eigenvalue of the corresponding total momentum operator, which commutes with the hamiltonian H. It 
is clear that p is an invariant to the "motion" of an eigenvalue as c varies continuously from 00 to —00. For each 
solution set {kj} there is another set {k'j} that also solves ( |l7| ) and is related to the former by 



k' 



27T710 • 



(23) 



The transformation {kj} — > {k'j} amounts to shifting the total momentum by 67m . This means that any state can 
be mapped to another state in the central momentum strip — 37r < p < Sir, and vice versa, by such a transformation. 
Thus we shall limit ourselves to study only those states having total momentum in this strip, namely, 77 p = —1,0, 1. 

A convenient set of variables, especially when the three fc's are real, that uniquely define the state is p, 81 and 82, 
where <5i and 62 are, respectively, the relative momenta between particles 12 and 23, namely 



<5i = k 2 
82 = h 



k 2 . 



(24) 
(25) 

Note that, when real, the order assumed for the kj, Eq. (pp[) , implies that Sj > 0. The fc's are given in terms of these 
variables by 



\(p-28 x -8 2 ) 



&2 = -^{p + Si - S 2 ) 



282 



(26) 
(27) 
(28) 



and the energy takes the form 



E 



-\p 2 + 2{8l+8l+8 l 82)} 



Combinations of Eqs. ( |l8| ) and (|l^) give the two coupled equations for 8\ and ^2, 



8\ = i In 
82 = i In 



+ \ c — iS 2 c + {(81 + 82) 



c — iSiJ + 182 — 1(81+82) 
c + iS2\ 2 c — iSi + 1(81+82) 



182) c + i8i — 1(81 + 82) 



27T771 



27T772 



(29) 

(30) 
(31) 
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where, on the basis that the principal part of the logarithm is taken, n\ is not necessarily equal to 777,2 — mi, nor is 
ri2 necessarily equal to 7773 — mi- In fact, unlike the rrij, the rij do not have to remain constant as a given "root" of 
the coupled equations {61(c), 62(c)} changes continuously with a variation of c. That is why we shall not classify the 
roots according to "local" values {ni(c), n2(c)}, but according to their values {n^n"} for no interaction, namely for 
c = 0. These are unambiguously related to the set of quantum numbers n%j, see Eqs. ( |35| ) and ( [36] ) below. If, as c 
changes, the argument of one of the logarithms, say Zj, crosses the negative real axis (which is the branch cut for the 
principal part of the logarithm), its phase changes abruptly by ±27r and the corresponding rij has to jump up or down 
by one unit in order to follow the root continuously. Of course, these discontinuities have no physical consequence and 
merely reflect the choice made for the branch of the logarithm. For several formal manipulations and in particular for 
the object of obtaining asymptotic expressions, it is useful to avoid the discontinuities by continuing analytically the 
logarithm across the branch cut, i.e., by passing to the contiguous Riemann sheet when the argument Zj crosses the 
real, negative axis. This will be discussed further in Sec. IV to clarify the trajectories of the kj as c changes smoothly. 

There are certain symmetries of the parameterization 61,62 and p which can lead to energy degeneracies. The 
interchange 61 # 62 together with the change in sign of p are equivalent to the changes fci # —k^ and &2 — fc'2 
which, for real fc's, inverts the order of Eq. (|2CJ) to —ki > — fc2 > —ka- This is also equivalent to taking the complex 
conjugate of the wavefunction. Thus, if k% ^ —k% and/or k2 7^ 0, these are two different states with the same energy, 
a two-fold degeneracy. On the other hand, if ki — —k 3 and /C2 = 0, which is the special case that 6\ — 62 and p = 0, 
then this symmetry reproduces the same state, the wavefunction is real, and the state is nondegenerate. In terms 
of the classification of the states of the central momentum strip by the 77°, it follows that interchanging n® and n 2 , 
(nf ^ TI2), amounts to the interchange of <5i and 62-, and, see Eq. (|37]) below, to change the sign of p, so that all signs 
of the fc's are changed and the complex conjugate state is obtained. But for the "diagonal" case, n° = nEJ, there is no 
degeneracy. These states are real and even under the parity transformation Xj — > —Xj. 

A second symmetry of the parameterization 6\, 62, p is the interchange 6± ^ — 62 while p remains unchanged. This 
is equivalent to the changes k\ ^ fc 3 and fc 2 k 2 , which inverts the order of Eq. (|2^) to k\ > fc 2 > k 3 but does not 
change the momenta themselves. But since the order of the fc's is immaterial, this is just another way of labeling the 
same state. In terms of the n° this means that (n\, n°) and (— n^, —n®) are actually the same state. 

Because of the stated symmetry relations and the fact that any eigenstate can be translated by the total momentum 
shift (E3|) to the central momentum strip, an exhaustive analysis of all possible states is achieved by examining the 
cases rv[> n\ > 0, since any other case is either equivalent to one of them or obtained by a simple transformation. An 
understanding of the behaviour of the roots and their limiting properties for different ranges of c requires a detailed 
analysis of how to carry out the analytical continuation of the logarithms in Eqs. ([}0|,|l]) as c varies. This is provided 
in the following sections. To keep track of the global picture a handy summary of the results is provided in the final 
Section, and a set of figures illustrate the essential aspects. 



III. THE REFERENCE CASE OF "NO INTERACTION" 

For c = there is no interaction and the particles move freely. In this case 

e lf ^ = -1 , (32) 
and with the present choice for the range of the Ojg, 

021 = 32 = 031 = -7T. (33) 

It follows that 

a(ijk) = 1 (34) 

5 1 = 27m? = 27r(m 2 - m x - 1) (35) 

52 = 2nn2 = 27r(m3 — 777,2 — 1) • (36) 



Equivalently, the kj are all real multiples of 2tt. Due to the conventional order (|20|), only the case where n° > needs 
to be considered in order to account for all the states of the system. 

The total momentum in the central momentum strip for a state (n®, n®) is determined by noticing, from (p7|), that if 
&2 is to be a multiple of 2tt, then p+6% — 62 has to be a multiple of 671. For any pair n\, n®, and p = 2nn v (n p = 0, ±1), 
there is only one possible solution for n p , namely 

if n\ — n® = in 

-1 if n\ - n\ = 3n + 1 n = 0, ±1, ±2, . . . . (37) 

1 if n\ - n\ = 3n + 2 
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For the states where one of the n° is zero, two fc's are equal. The equality of two fc's can occur only at c = and 
at certain critical c values discussed in Sec. V. (In general, for c =/= 0, the wave function vanishes if two kj are equal.) 
For the ground state, n\ = n 2 — 0, the three kj are equal and the wave function is a constant. 

The classification scheme used in the remainder of this paper, for the energy eigenvalues and states, is based on 
following Si and S 2 as continuous functions of c from the reference non- interacting system. A symbol such as (1,2) 
gives the values of the quantum numbers and n 2 and identifies a given "root" {<5i(c), S 2 (c)} of the transcendental 
equations and the corresponding eigenstate (within the central strip of total momentum) irrespective of the value of 
c. Note that the three "quantum numbers" equivalent to the set of quantum numbers {rrij}. The total 

momentum given in (^) is independent of the value of c, as well, the degeneracy of a root is also an invariant to the 
"motion" of the root with c. 



IV. REAL k SOLUTIONS 

On the assumption that the kj all remain finite (and real), it follows that 

if c — > 00; 

— 2tt if c — > —00, 

V31 > 

with the consequence, see Eq. (|l8|), that 




(38) 



— 00 < — c — ► 00 
27r(mi + 2) <— k\— > 2i:mi 

2irm 2 <— Al2^ 2lXm 2 

27r(m 3 - 2) <- fc 3 -» 27rm 3 . (39) 

Bounds on the A; differences can be narrowed by examining separately, the detailed properties of the (9's for positive 
and negative c, by the following reasoning: On the basis of the chosen order for the kj, Eq. (p(il), it follows that 

k 3 - ki > k 2 - h, fc 3 - k 2 . (40) 

As a consequence, for c > and — n < 9jg < for (j > £), 

««»(%)<««»(>),<.«..(>) (41) 

3 i < 21l 6 32 , (42) 
- 7T < 6 31 - 9 32 < (43) 



and 



so that 



and 

27rrii — 7r < 2Tr(m 2 — mi) — n + 29 2 i < 81 = k 2 — k\ = 2i:{m 2 — mi) + 29 2 i + 6 3 i — 9 32 

< 2-n{m 2 - mi) + 26 2 i < 2i:(n° + 1). (44) 

An analogous argument for c < leads to 

27r(n? -l)<Si = k 2 -ki< 2itn\ + tt, (45) 



with the lower bound approached according to Eq. (g9(). Bounds for S 2 involving n 2 have the same structure. 

It is thus seen that if n® > 2 and n 2 > 2, then these k differences remain positive as c varies from 00 to —00, with 
the consequence that the kj are real for all c under these conditions. This also implies that the energy tends to a 
constant value for \c\ — > 00, see Fig. 1. For these states, the delta function potential acts, for both very large positive 
and negative c, effectively as an impenetrable barrier. That is, the wave function at the region boundaries, Xi — Xj, 
tends to zero. This can also be deduced from the jump conditions, Eq. M), since the derivatives produce only the 



G 



finite kj while c becomes infinite, so that consistency requires that ip vanishes. The potential energy as a function of 
c decreases from zero to an intermediate minimum and then grows again towards zero as c — > — oo and vice versa for 
c > 0, see the Appendix and Fig. 2. 

The above deduced limitations on how the 5's change with c allow a more detailed discussion of these changes from 
the point of Eqs. ( f30| ) and (31). An appropriate starting point is to consider the limit of these equations when c — » 0, 
namely 



5% = i In 



i In 



8\ — ic\ &2 + ic (5\ + 82) — ic 



Si 
62 



82 — ic (81 + 82) + ic 



ic\ 81 + ic (Si + S 2 ) 



ic 



81 — ic (Si + 82) + ic 



27771? 



27T7J 



2 ' 



(46) 
(47) 



where the factors in the arguments have been rewritten to emphasize that the <5's are in this limit the leading terms. 



The identification of the nj is on the basis that, on expanding these equations as c — > 0, Eqs. ( p5|) and ( p6|) are 
obtained, namely 



61 = 2im\ + 



2ngn| + 2(n°) 2 - (n\f 
n\n f 2(n\ + n^ir 



(48) 



with a symmetrical expression for 82 on interchanging and n°. These expansions are valid only if both nP- 7^ 0, 
while the exceptional cases are examined in Sees. VI and VII. 

As c increases positively [always assuming the order of Eq. (pfj)], the phases of the factors in the arguments of 
the logarithms change, compare the Bj t of Eq. (H) and the discussion of their behaviour. On tracing this behaviour 
as c — > +00, it is seen that the phase of the argument _of each logarithm decreases by 27r, so this can be taken into 

2 + 1. For negative c 



account when changing Eqs. ( |46[) and ([47J) into Eqs. (BOh and fl3l|) by setting rij(c — > cxd) 
the phase changes are the opposite and Eqs. (Eq) and (MTj) are appropriately changed into 



Si = i In 
82 = i In 



i(5i\ — c + i<5 2 — c — + $2) 
-c — %82 -c + 1(61+62) 



iSi 



—c — 182 \ " —c + i8\ —c — i(5i + 82) 
-c — iSi —c + i(Si + S 2 ) 



iS? 



+ 2w(n° 1 - 1) 
+ 2n(nl - 1) 



(49) 
(50) 



These are appropriate for expansions when c —> —00, but may of course be used for all c by analytic continuation. In 
the same vein, the related pair of equations with nj = rz5 + 1 may be regarded as valid for all c by analytic continuation 
across the logarithm branch cut. 

In summary, provided both n® > 0, then 



+ 1) (l - -\ +()(,-). 



(51) 



while, provided both > 1, then 



3 27r(n°-l) I 1 + — ] + 0(c~ 2 ). 



(52) 



In contrast, if n\ = 1, then 81 can approach for finite negative c, see Fig. 3. For more negative values of c, the kj's 
can become complex. Similarly for = 1. Section V discusses this situation. If one of the n° = 0, then 8j vanishes 
at c = 0, see Sec. VI. 

The case in which n\ = is particularly simple to analyze. It follows from Eqs. (flfj]) and (^), analytically 
continued for all c, that Si and 82 satisfy the same equation and thus are equal, with the consequence that the three 
k's for such an eigenstate remain equally spaced as c varies. After dropping the subscripts, and formally written for 
c > 0, the equation for the common 8 is 



8 = i In 



(c + iS)(c + 2iS) 
(c~iS)(c-2i5) 



+ 27r(n° + l). 



(53) 
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Fig. 4 illustrates how 8 varies with c. By implicit differentiation, 

dS = 6<5(c 2 + 2<5 2 ) 

dc c 2 (c 2 + 5<5 2 )+4<5 4 + 6c(2 ( 5 2 +c 2 ) ' 1 ' 

If n\ = n\ ^ 0, Eq. (^||) is consistent with the c — > limit of this result. 

V. STATES WITH AT LEAST ONE n° = 1; DIMER STATES 

If one rij = 1, for definiteness nf = 1, then from Eqs. and (fl5|), there is a possibility that <5i becomes for 
some critical negative value (7(1, 71°) of c. Note that the equations for 8\ and 82 are independent of p, so the critical 
value C(l,7i2) has nothing to do with the value of p. [The discussion for the degenerate partner (n°, 1) follows similar 
lines substituting n\ by nS], <5i by 82 and p by — p, see the discussion of the symmetries of this parameterization in Sec. 
II.] At a critical point (7(1, n®), fci and k 2 become equal (if = 1 as well, all three k's become equal) and the Bcthc 
ansatz (|l^) form for the wavefunction is no longer valid because it vanishes. Of course the normalization constant also 
vanishes so that the normalized wavefunction does not vanish, but merely has a different functional form obtained by 
taking the limit as c — > (7(1, n®) using the rule of l'Hospital. A similar case was found by Lieb and Liniger for one 
particular root in the case N = 2 || . 

Two different cases arise, according to whether n 2 > 1, or if n\ = n 2 = 1- These cases are discussed in turn. 



A. The case when n§ > 1. 

The value of the critical point (7(1, n 2 ) and the behaviour in its neighborhood can be obtained by expanding Eqs. 
(||) and © for 81 -> 0. For such a purpose it is useful to introduce the factors Uj = Sj/c, j = 1,2, and define 
uq = 112(0 = C) as the critical value of U2- After fairly extensive algebra, it follows that 

with 

C(l,n° 2 ) = -A--^, (56) 
1 "+■ u o 

where uq is determined by the transcendental equation 

3i lni±^+27r(n°-l)+4u + -^=0. (57) 

1 — IUO 1 + Uq 

In the neighborhood of the critical point, U2 changes according to 

1 3 + 6m^ + 2u* 2 

u 2 = uo - ~u x + - 2 ° u\ + ■■■ . 58 

2 6mo(1 + Uq) 

It is seen that Si(c) has a square root singularity at the critical point, being real for c > (7(1, n§) and pure imaginary 
for c < (7(1, n§). An analytic connection between real and imaginary branches of S\ can be made by attributing c 
with a small imaginary part. The association that is used in the following parameterization is consistent with c having 
a small negative i mag inary part, essentially adding a small dissipative contribution to the Hamiltonian. It is also 
noticed from Eq. ([36|) that all the critical values of c lie between —6 and —4. Actually, to three decimal places, the 
lowest critical value is (7(1, 2) = —4.163, increasing towards —4 as Uq, and n®, increases. These aspects are illustrated 
in Fig. 3. 

For c < (7(1, n§), fci and &2 become a complex conjugate pair, while remains real. The discussion of the root in 
this region is better examined by using a new set of real variables, a and 7, defined by 

81 = —Ha 

82 = ict — 37 , (59) 
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so that 



ki — ia + 7 + p/3 
fc 2 = -ia + 7 + p/3 
fc 3 = -2 7 + p/3. 



(60) 



This is consistent with Eq. ( J5q ) for a — > 0, on the basis that both a and 7 are real. This parametrization could 
of course be used for all c [with a possibly imaginary], since 8\/2 = —ia has the physical meaning of the relative 
momentum between particles 1 and 2 while 27 is the relative momentum between the pair 12 and particle 3, 



5i/2 = 



k 9 - fci 



p 2 / fci + fc 2 
2 7 =fc3-3 = 3(fe3-^— 



(61) 
(62) 



For real a, l/a gives a measure of the size of the "bound states" formed. Of course, in a finite box all states are, 
strictly speaking, bound, i.e., their energies are discrete and their spatial extention is limited by the box length. But 
when a is real, the state is localized even further so that the probability of a pair of particles being close together 
has been significantly enhanced. As well, the energy has a negative contribution, see ( p3[) below, so that an energy 
gap arises in the spectrum between states with real or imaginary a, see Fig. 1. Since these are all basic ingredients 
of proper bound states, this terminology seems justified. 

With this parametrization, the total energy is decomposed into separate quadratic contributions from the three 
variables, 



E 



-2a 2 + 6 7 2 + p 2 /3, 



and the system of transcendental equations takes the form 



2a = In 



37 = i In 



2a 



-c — a — 3*7 —c — a + 3*7 



z + 2a, 
-c + a + 3*7 



c + a + 3*7 — c + a — 3*7 
2 



-c+ 2a 



3*7 



a — 3*7 / —c — 2a — c + a — 3*7 



+ 2 7 r(r* 2 l - 1) 



(63) 

(64) 
(65) 



which corresponds to the correct phase form for c < 0, according to Eqs. (49) and (pd) 



It is consistent to solve these 

equations maintaining a and 7 real, which is also consistent with the local behaviour, Eq. (58), as a — > 0. Before 
entering into the detailed analysis of the equations, it is worth examining Figures 5 and 6 to quickly visualize the 
behaviour of these two parameters with c. Solid lines correspond to n\ = and dashed lines to n\ — 1, whereas 
the numbers close to the different lines give n®. As c becomes more negative a increases. This concentrates the 
wavefunction to where the particles are close together and (V) and (E) become very large and negative, see in Figs. 
1 and 2 the lines for n — 0,1. A prominent feature in Fig. 5 is the grouping into two asymptotic behaviours for 
aasc-> —00. These will be later associated with dimer and trimer configurations. Note also that all states with 
n\ = have a common critical point at c = (where a = 0), while for rtj = 1 the critical points spread from c = —6 
to c = —4. Another interesting point is the quasi-invariance of a with respect to n\ for > 2 at fixed c. These 
states have essentially the same binding strength and differ only by pair-single relative momentum, and possibly by 
total momentum. Fig. 6 for 7 has a simpler structure based on a rather regular pattern of lines with equally spaced 
asymptotic values. Except for the two cases (0,0) and (1,1) with 7 = for all c, which corresponds to no relative 
motion between a particle pair and a single particle, when c becomes more negative, 7 varies smoothly and tends to 
a constant value. When 7 is essentially constant, c changes the strength of the attraction between the pair (i.e., the 
value of a), but not the motion of the third particle with respect to the pair. 

The structure of the spectrum of energy levels can be described in terms of how the energy varies as the three 
quantum numbers n\,n^,n p change. This can be attributed to several types of "elementary excitations" which are 
associated with different physical effects. Within the central momentum strip, for states with complex kj's and for 
— c large, these are: 

a) As r*2 varies, a change of relative pair-single motion by A 7 w 7r/3, with a essentially constant and \p\ constant. In 
Fig. 6 this is not possible between all contiguous levels of 7, but only for those where n p = 1 and n p i = — 1, see Eq. 
(37). If one of the states has momentum zero, a jump to the nearest level necessarily implies in addition to A 7 , an 
elementary total momentum jump A p = 27r; b) Transitions between a trimer and a dimer state or from a dimer to a 
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pair-absent state, with 8 a ~ — c/2; c) A minimum total momentum jump by 2tt, with a and 7 constant. This may 
only occur between the states (0,0) and (0,1) or (1,0). Any other transition changes 7. But a and 7 may also stay 
constant il the system changes to a different momentum strip by a total momentum translation of A p = 6ir. Of course 
multiples or combinations of these elementary excitations are possible and complicate the spectrum considerably. 

The detailed quantitative features of a and 7 as functions of c are now examined. On the basis that a and 7 are 
real, the real part of Eq. ( |65|) is 



3., 
—i in 

2 



(— c + a — 3i7)(— c — a — 317) 
+ 3i 7 ) 



(— c + a + 3i7)(— c - 
This provides upper and lower bounds for 7, namely 

- (4n° - l)7r/6 < 7 < -{An% - 7)tt/6 

Equation ( |66|) may also be written as 

67c 



= 37 + 2tt(^ - 1). 



arctan 



9 7 2 



-7 - 2tt(^ - l)/3 



(66) 



(67) 



(68) 



which is useful for the determination of the asymptotic behaviour of a and 7. 
It is also possible to solve for 7 2 in Eq. (H) in terms of a and c, 



^ 2 _ er a (c-2a)\c - a) 



e a {c + 2a) 2 {c + af 



9[e«(c + 2a) 2 -e-«(c-2a) 2 ] 



(69) 



As a — » 0, 7 2 (a = 0) = c 2 (6 + c)/[— 9(4 + c)]. This is identical to Eq. (56) for the relation between 82 and c at a 



critical value of c. 

The asymptotic behaviour as c — > —00 is now investigated. According to Eq. (p^), 7 remains finite while a 
satisfies Eq. ([64]). Both positive and negative a are solutions to this equation, but since these just correspond to an 
interchange of ki and hi, only the positive root is examined. Eq. (64) can be rewritten as an equation involving only 
real quantities, namely 



2a = In 



-c-2a\ (-c- a) 2 +9j 2 



-c+2a/ (-c + a) 2 + 97 2 



(70) 



The obvious (but invalid) approach to try when making an asymptotic expansion, is to expand in powers of 7 2 since 
the other factors involve c. This implies that a must also remain finite and leads to the requirement that c — ► —6, an 
inconsistency. It follows that either (— c — 2a) — > or (— c — a) — > in the Limit c — ► — 00. But it is noticed that as c 
changes from its critical value Cl\,n^) to —00, a changes from to its asymptotic behaviour, yet the factors in the 
argument of the logarithm, — c ± 2a and (— c ± a) 2 + 9"f 2 , must remain positive or an imaginary phase factor must be 
added to the right hand side. Since such a case would imply that a becomes complex, this is not allowed. Another 
way to understand the preservation of sign of all factors is that if one of them became zero for a finite c (and a) the 
logarithm, and the left hand side of the equation, would be infinite in absolute value, which is again inconsistent with 



the finite value of a on the right hand side. The only form of a that maintains all factors positive, 
with (3 < and (3/c — > asymptotically. A straightforward expansion then gives 

(3 = 3ce c/2 - 9c 2 e c + • • • , 



a 



-c + 'ice c/2 - 9c 2 e c + ■■■, 



■0, 



(71) 



and from Eq. (pq). 



7 = ~(«8-l) 



1 - - 

c 



(72) 



as the asymptotic expansions for both a and 7. 
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B. The case when n° = 71? = 1. 



As previously mentioned, this implies that Si = 82 for all c, and the common 8 is determined by Eq. (|5 
appropriate for c > 0. Reexpressing this for c < 0, retaining a continuous relation for the phase, gives 



i In 



(-c-iS)(-c-2i8) 



{-c + iS)(-c + 2i5) 



which is 



(73) 



and taking into account that n° — 1. Clearly <5 = is a solution of this equation. But to find the corresponding 
critical value C(l, 1) of c and the behaviour in the neighborhood of this critical point, an expansion is needed. This 
is easily accomplished to yield 



-6+^ 2 + - 
6 



(74) 



This immediately shows that the critical value is C(l, 1) = —6, and that there is a square root singularity of S as 
a function of c. For c < —6, 5 = —ia is pure imaginary [negative imaginary if the same connection around the 
singularity is used as in the last subsection] . It follows that the kj are given in terms of this parameterization by 



&3 = — ia + p/3 

k 2 = p/3 

ki = ia + p/3, 



while a is determined by 



In 



1 - 3a' + 2a n 
1 + 3a' + 2a'' 



(75) 



(76) 



where a' = a/c. [Contrast this with Eq. (60). The associations made for the kj in (£51) are more natural here, since 
Si = 82 and the three kj change continuously across the critical point, but in fact (|60D, with 7 = could be used as 
well because the state is defined by the three "momenta" regardless of the ordering convention.] 

Clearly if a' — > as c — > —00, then a — > 0. But on expanding to look at the correction terms, this assumption 
implies that c — > 6, an inconsistent result. It follows that a grows to 00 asymptotically. On the basis that only 
a > solutions are needed (a' < 0), the vanishing of the denominator in the argument of the logarithm requires 
either a' — ► —1 or a' — > —i as c — > —00. Only the latter is consistent with no accumulation of phase as c changes 
from —6 to —00, and the reality of a. (Note that, by a similar argument to the one below Eq. J7C|), numerator and 
denominator in (|76|) must preserve their sign as c varies from —6 to —00 which implies the bound 2a < — c.) On 



setting a 



7], it follows on expansion that 



a = --c + 3ce c/2 



(77) 



as the asymptotic behaviour of a in this case. 



VI. STATES WITH ONE n) = 0: DIMER AND TRIMER STATES 

This group of states are of the form (0, nSJ), and of course their degenerate partners. Only the case (0, nj]) is 
explicitly treated here since their corresponding partner states are simply obtained by symmetry, specifically for the 
states when c < by changing the sign of 7. For c — > 00, these states fit into the formulation of Eq. (pl|), so no 
further discussion of this limit is needed. 

As c — > 0, the starting point for the analysis is the pair of equations, ( |46| ) and (f47j). Since n\ — 0, 8\ and c are 
to simultaneously approach 0, and it is to be expected by analogy with the other cases, that c will be proportional 
to 8\. This is confirmed by the following argument: Since 8\ — > 0, then the argument of the logarithm in Eq. ( |46| ) 
must approach 1. For the ratios of the terms involving 82, these approach 1 since they are dominated by the common 
non-zero value of 82- That leaves the square term. For this to approach 1, each factor must be dominated by the 
common Si, which implies that it is the ratio c/Si that is small and can be used as a variable for which the logarithm 
is expanded, as well as the ratios 0/82 and 8\j '82, all of which vanish as c — > 0. Explicitly, the expansions of Eqs. ( |4^ ) 
and @ are 
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Si 



4C3 2^ + 0(dt) 



5i 35? 5 2 2 



c 2<f 
% + 35? 



5-. = 2rr„!J - 2- ^+6^-2^ + 0$) 



(78) 



It follows that c is proportional to Sf, as was to be deduced. Rearrangement of these series for c — > gives the 
expansion in powers of 5i as 



1 



1 



192 32(7m°) 2 



r o Ir 35? 



(79) 



For c < 0, the parameterization of Eq. ( |59| ) is appropriate. In the limit c — > —0, the ratio c/8\ transforms according 
to 



c 



1 



1 c 
-—ia <— 17—, 

2 2a' 



with the consequence that Eqs. (Efih and ( p7| ) become 



2a = In 



ia — 37 — i In 



2a + c\ (a + c) 2 + 97 2 



2a- cj (a - c) 2 + 97 2 

2a — c\ f ia — ic — 2>^\ l f —ia — ic — 3-f 
2a + c J \ia + ic — 3j J \—ia + ic — 3'j 



2nni 



(80) 

(81) 
(82) 



It is also important to note that, as a consequence of the analytic continuation into the c < region, 2a + c > 0, 
and 7 — > —^2/3 < 0. To avoid further singularities in the logarithmic expressions (which would lead to inconsistency 
between right and left sides of the transcendental equations), these constraints must hold for all c < 0. 

For exploring the behaviour as c — » -co, Eq. (pl| ) for a has an appropriate form for expansion about the large 
—c and a. But if one term in each of the factors is to dominate the expansion, then the result would imply that the 
logarithmic expression becomes finite, a result inconsistent with the left hand side approaching —00. As a consequence, 
since a > — c/2, either a — > —c/2 or the combination, a — ► — c and 7—^0, must occur. The first alternative is now 
shown to be valid only for n® > 1, while the second alternative is valid only if n® = 1. 

The limiting case a = — \c + (3 leads to a straightforward expansion of Eq. (jsi|), 



='--' = -"' I 4:) + l£ 



(83) 



which can be rewritten as the equation 



/3 = — 3ce 



c/2 _ 



9cV 



(84) 



for /3. As c — ► —00, this vanishes exponentially. The behaviour of 7 is to be obtained from Eq. (|82j). But the dominant 
quantity c is multiplied by the phase factors ±i, so that an asymptotic expansion carries along a phase change for 
the logarithm. On carefully analyzing how the various factors change as c changes from to -co, a phase change of 
e 37 ™ in the argument of the logarithm is found. The real part of Eq. ( |S2| ) determines 7, which after expanding and 
rewriting gives 



" = - ( 3^ - 1 



1 

c 



(85) 



Since it is required that 7 < 0, it is seen that this expansion is only valid for > 1. 

The limiting case a = — c + r\ also allows a straightforward expansion of Eq. (^l]), but in this case the resulting 
equation for rj involves 7 in the lowest order term, namely 



r] 2 + 7 2 = 36e~ 2,? c 2 e 



2c 



5r/ 
3c 



(86) 
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Since both 77 and 7 are real, this implies that both these quantities must vanish asymptotically as c — > —00. In this 
case, Eq. (|82| ) accumulates a phase e 37 ™/ 2 in the argument of the logarithm when transforming the factors so that 
they will be dominated by a positive real part in the limit. What is crucially different in this case from the previous 



one is the presence of finite complex factors —37 =p ir\. As a result, the real part of Eq. (82) has the asymptotic 
expansion 

|axctan(^) + (2nO-|) 7 r-| 



3 7 = - arctan [ ^ ) + I 2n° 2 - - ) n - + ■ ■ ■ . (87) 



As c — > —00, 7 must vanish for this case, as was deduced from the expansion of the a equation. Thus in the limit, 
the identity 

V \ _ f 4 „ 



lim arctan = -n, — 1 

must be satisfied. Since the magnitude of the arctangent is bounded by 7r/2, this identity can only be satisfied if 
n® = 1, in which case the limiting ratio of 77 and 7 is determined by the condition 77 = — 37 tan(7r/3). In summary, 
for nSj > 1 the asymptotic behaviour of a is a — > — c/2, which gives states of dimer type, while for n\ = 1, a — > — c 
and all three particles are forced to be close to one another, a trimer state. These associations will be discussed in 
Section VIII. 



VII. THE GROUND STATE: n? 



n 2 



Since and n° are equal, it follows that 8% = 82 = 5- Near c = and for c > 0, the appropriate equation for 
determining 8 is Eq. (Eq), which is modified for equal S's to be 



i In 



(<J-ic)(2<J-ic) 



(5 + ic)(2<5 + ic) 



(89) 



The behaviour as c — *• 00 is covered by the expansion of Eq. (p3|), with asymptotic form identical to Eq. (|5l|). Since 
5 — > as c — * 0, it is necessary that the argument of the logarithm must approach 1, which requires that c — > faster 
than does 8. Thus the expansion parameter is c/8, so that after expansion and rearrangement 

C = ^ + I L S « + .., (90) 

For c < 0, the parameterization of Eq. ( |75| ) is used, this being equivalent to 8 = —ia with a > 0. The ratio —ic/S for 
c > thus becomes c/q for c < and Eq. ( |89| ) becomes 



= In 



(a — c)(2a — c) 



(a + c)(2a + c) 



(91) 



This is identical to Eq. (|76|), but now the constraint is that a > — c. This constraint requires that as c — > —00, a 
must become infinite and the expansion of the logarithm is about a singular point of the logarithm. The only possible 
form is a = — c + 77, with 77 > approaching zero. After rearrangement, the resulting expansion gives 

?7 = -6ce c - 36c 2 e 2c + • • • . (92) 

which is consistent with 77 being positive. 

VIII. TYPES OF STATES AND THEIR REPRESENTATION 

The states are best represented as contour plots of the probability density (and of the phase if required) in a 
"ternary phase diagram" for the variables 

ri2 = x 2 - x\ 
T23 = x s - x 2 

rsi = 1 + Xi - x 3 , (93) 
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constrained to the region -R123, < ry < 1. Note that these three coordinates always add to one. In this diagram 
the coordinate points are represented in an equilateral triangle. Each of the base lines corresponds to one of the 
coordinates being zero, and each point in the base line corresponds to a particular location of the third particle 
at the right or left side of the pair ij (the closer to the vertex, the closer the third particle is to the pair). The lines 
parallel to the base are lines of constant . The value of r*y increases from zero at the base to one at the opposite 
vertex (labeled as in the figures). The center of the triangle is the point where the three distances are equal to 
1/3. Near the vertex r%j the distance between particles j and i is also small (and tends to zero at the vertex itself); 
the difference with the basis region is that now the third particle is between the particles j and i. 

In summary, bases are associated with two particles being together (dimer configurations), and vertices with the 
three particles being together (trimers). However, in a general state with three real k's, there is no bias towards 
these configurations. Recall that the wave functions in -R123 are linear combinations of six exponentials that can be 
obtained from the (123) form, 

e i(kix 1 +k 2 X2+k 3 x 3 ) (Q'i) 

by permuting the fc's in all possible manners. The probability density (square modulus) of any of these real k plane 
wave terms is constant; in other words, in these plane waves none of the particle configurations is favored. The 
interference between the six different plane waves however destroys the spatial homogeneity (except for the ground 
state at c = 0) and provides some structure with maxima and minima, see Fig. 7. The complex k case is different, 
see Figs. 8 and 9 for examples of trimer and dimer states. Eq. (|)4|) for the (123) exponential may now be written, 
using the parameterization in (|60|), as 

e i(x 1 +x 2 +x 3 )p/3 e ii(x 1 +X2-2x 3 ) e a(x 2 ~x 1 ) {95) 

where the plane wave for the center of mass motion, a plane wave for relative motion of the pair 12 with respect 
to particle 3, and a real exponential can be recognized; for a > the exponent is positive and it favors the trimer 
configuration of the vertex (r 12 — 1). By permuting ki and k 2 , one finds instead, in the term (213), a negative exponent 
that favors the dimer configurations ri2 = 0. Of course the other two pairs, 13 and 23, have also a corresponding 
set of dimer and trimer contributions, so that the six terms of the wave function can be separated into two groups: 
Three of them, (123), (231) and (312), represent trimer configurations, and the other three, (213), (132) and (321), 
represent dimer configurations. The relative weights among them are determined by the amplitudes a^fc. For "trimer 
states" the three trimer terms dominate the linear combination and the energy becomes, as c — > — oo, the energy of 
an actual trimer state (for three particles on an infinite line). For "dimer states", there is also significant density 
along the edges of the triangle (not only at the vertices), and the energy tends to the energy of the actual dimer (on 
the infinite line) , plus the contributions from relative motion of the dimer with the free particle and of the center of 
mass motion. This is consistent with our expectation of reproducing infinite line results in the limit of a large box. 
As a concrete example, the (unnormalized) states (0,0) and (1,1) for c < —6 are examined: In both cases 7 = p = 
and they can be written, using Eqs. ([l^), (|T4|), ( |ilf ) and (pl|), as 



if, = e ari2 + e 
2a - c 
2a + c 



(e~ ari2 + e- Qr23 + e- ar31 ) . (96) 



As c — > —00 the factor multiplying the dimer terms in parenthesis tends to 3 for (0, 0) (which makes this contribution 
negligible) but to 00 for (1, 1). Fig. 8 shows the ground state (0,0) for negative c. 

An important aspect of these associations is that the trimer or dimer character changes continuously along a given 
root as c varies, and only asymptotically (c — > —00) is the separation between trimer and dimer states unambiguous. 
For any finite negative c the complex k roots have non-zero dimer and trimer components. Note for example how the 
state (0,1) goes from a dimer dominated behaviour to trimer behaviour as c becomes more negative in Fig. 5. In the 
same vein, even though the threshold between real and complex k is well defined and it occurs at a critical value of 
c, there aren't any dramatic (discontinuous) changes in the wave function, and the energy varies smoothly with c, see 
Fig. 1, in the neighborhood of the critical c values. However, a different qualitative behaviour (of the energy and state 
probability density) becomes clear when comparing the state below and above the critical point as the distance from 
C increases. Thus the critical values indicate a transition of the root from one character, without pair formation, to 
another where dimers or trimers can be recognized. 



14 



IX. SUMMARY 



A model of three bosons subject to delta function interactions and periodic boundary conditions has been analyzed. 
In particular a description of the eigenstates and their behaviour has been given in terms of three momenta fcj, 
i = 1,2,3 or two sets of alternative parameterizations, Si,d2,p and a,j,p, convenient, respectively, for the cases 
where the fc's are real or complex. The roots can be primarily classified according to whether the three momenta 
remain real, or not, for all c. In the second case the wave function tends to concentrate asymptotically around dimer 
or trimer configurations and the energy decreases quadratically with c as c — > — oo. The critical values of c required 
to form the bounds (go from real to complex fc's) have been provided. 

The main features of the root behaviour as c varies are now summarized. (Using the symmetry properties or total 
momentum translations the behaviour of any other state is obtained from the ones we consider explicitly, namely 
states (nl, n®), n± < n^, in the central momentum strip, n p = 0, ±1.) Relations satisfied by all roots are: 



The total momentum p (within the central strip) is given from n\ and n® by (|37|) . It is constant as c varies for 



a given root. The total energy varies according to (29) or (p3J). 



• Sj{c = 0) = 27m°. 

• Sj — > 27r(n^ + 1) as c — > oo. 

The different particular cases are characterized by the following properties: 

• (0,0): C(0,0) = 0, p — 0, 7 = 0, a ~ — c as c — * — oo. The ground state is a nondegenerate state with trimer 
character. 

• (0,1): C(0, 1) = 0, a ~ — c and 7 — > as c — > —00. Asymptotic trimer character. Similar to the ground state 
but it has a degenerate partner and p 7^ 0. 

• (0, n\ > 1): C(0, 7i§) = 0,7^ (— 2/3n,2 + 1)tt and a ~ — c/2 as c — > —00. Dimer character. 

• (1, «■? > 1) : ~6 < C(l, n°) < —4, 7 — > — 27r/3(ri2 — 1), a ~ —c/2 as c — > —00. Also dimer character, but it takes 
a stronger interaction to achieve in comparison to the previous group. 

• (n\ > Trirj — n i) : R- ea l ki for all c, <5j(±oo) = 27m° ± 1. The energy tends to a constant value as \c\ — > cxd. 

It is hoped that the root structure found and the techniques developed for its study will be useful for the examination 
of variants of the model involving different interactions and/or an arbitrary number of particles. Since all eigenstates 
can be obtained for a given c, one of the possible applications of the present work is the simulation of time dependent 
rearrangement processes using a discretized basis. In this context it may serve as a reference exact model to compare 
with approximate wave function propagation methods based on periodic boundary conditions 0. The model may 
also permit an explicit comparison between the classical concepts of "bound pair" and "free" particles and their 
collisional rearrangements with their quantum counterparts. The origin of the difficulties in the quantum case is that 
the Hamiltonian for the 3 (or N) body system does not commute with the Hamiltonian of a pair, so that using the 
concept of the bound pair in a gas (or in a box) is a delicate matter |3l| ]. 
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APPENDIX A: NORMALIZATION AND POTENTIAL ENERGY 

The expression ( |l2"| ) for the wave function is not normalized. For the analysis of the energy spectrum this does 
not have any effect. However the calculation of physical expectation values, and in particular of the potential energy, 
requires a normalized function. This requires the evaluation of the inner product (ip\ip), where \tp) is the unnormalized 



total wave function given explicitly by Eq. (12) in the region i?i23- The contribution of the six regions is identical so 
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(ip\ip}=6 / dx 3 / dx 2 / <ia;i|^(fci,fc2,A;3;a;i,a;2,a;3)| 2 (Al) 
Jo Jo Jo 

This integral is decomposed into 36 terms with integrals of the general form 

/ dx 3 / da; 2 / dx ie t{aiXl+a2Xl+a3Xl) (A2) 
Jo Jo Jo 

where the en's are combinations of momenta of the form — k* + ki, = 1,2,3). These integrals can of course be 
solved explicitly but the result is so lengthy that it is not reported here. In order to handle all terms efficiently, it is 
useful to classify the possible types of integrals. There are three cases: (ctj =0, j = 1, 2, 3), (aj — 0, ak + a.i = 0), 
and (oij ^ 0, j ' = 1, 2, 3). In all cases the sum of the en's is zero, J2j a j = P — P = 0. 

Once {ip\ip) is obtained, the (dimensionlcss) potential energy involves a simpler integration because of the delta 
functions, namely 

6c f 1 f X3 

{V) = Jo dXz J dx ^( x ^ x ^ x s)\ - (A3) 

Again, all integrals involved can be carried out. 
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FIGURE CAPTIONS 

Fig. 1 (E) vs c for (n ,n ). From bottom to top, n° — 0, 1,2,3. The lines for n° — and 1 extend to — oo as 
c — > — oo. 

Fig. 2 (V) vs c for From bottom to top, n° = 0,1,2,3. The lines for n° = and 1 extend to — oo as 

c — > — oo. 

Fig. 3 Si vs c. The numbers close to each line correspond to n® and n®. The figure for 62 is identical by 
interchanging n\ and n!]. The lines are only drawn up to the critical values of c where 8\ = [for (1,1), (1,2) and 
(1,3)], or 82 = [for (2,1) and (3,1)]. For more negative values of c, the a-7 parameterization is used, see Figs. 5 and 
6. 

Fig. 4 5 vs c for (n°,n°). From bottom to top, n° = 0, 1, 2, 3. 

Fig. 5 a vs c. Solid lines: = 0; Dashed lines: nj — 1. The number close to the lines is n°. 
Fig. 6 7 vs c. Solid lines: n° = 0; Dashed lines: = 1. The number close to the lines is n°. 
Fig. 7 Contour plot of the probability density of the state (3,3) at c = 1. The interpretation of a "ternary phase" 
type of diagram is explained in the text. 

Fig. 8 Contour plot of the probability density of the state (0,0) at c = —9. 
Fig. 9 Contour plot for the probability density of the state (0,2) at c = —9. 
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